arXiv:1503.08753vl [nucl-th] 30 Mar 2015 


Compact Stars in the QCD Phase Diagram IV (CSQCD IV) 
September 26-30, 2014, Prerow, Germany 

http: //www. ift .uni . wroc .pl/'^csqcdiv 


Core-collapse supernova matter: light clusters, 
pasta phase and phase transitions 


Helena Pais^, Constanga Provideneia}, William G. Newton^ and Jirina R. Ston^'^ 

^ Centro de Fisiea Computaeional, Department of Physics, University of Coimbra 
P-3004-516 Coimbra - Portugal 

“^Department of Physics and Astronomy, Texas A&M University, Commerce, TX 
75429-3011, USA 

^Department of Physics, University of Oxford, Oxford 0X1 3PU, United Kingdom 

"^Department of Physics and Astronomy, University of Tennessee, Knoxville, TN 
37996, USA 

1 Introduction 

The complex structure of nuclear matter in the density region approaching pg ~ 0.16 
fm“^ (central density of heavy nuclei) at finite temperature (T < 20 MeV) critically 
affects many astrophysical and nuclear physics phenomena. At low densities, the 
frustrated system called the “pasta phase”, caused by the competition between the 
Coulomb interaction and the strong force, appears, and is constituted by different 
geometrical configurations, as the density increases pQ El El IH El El El • 

The main interest in the pasta phase in core-collapse supernovae (CCSN) is that 
the neutrino opacity, which plays the main role in the development of a shock wave 
during the supernova collapse, is affected by its presence |H]. At very low densities 
(up to 0.001 times the saturation density), light nuclei (deuterons, tritons, helions, 
«-particles) can appear, and, like the pasta phase, can modify the neutrino transport, 
which will have consequences in the cooling of the proto-neutron star [HI EDI ED E2] • 
In this work, we use different methods to study the pasta phases of CCSN matter: 
a non-relativistic method, the 3DHFEOS method, where a 3D Skyrme-Hartree-Fock 
approximation is used, and relativistic mean field (RMF) methods, the Thomas- 
Fermi (TF) and the Coexisting Phases (CP) approximations, and a thermodynamical 
spinodal (TS) calculation. 

2 The 3DHFEOS method 

We use a 3D Hartree-Fock approximation with a phenomenological Skyrme model 
(3DHFEOS) for the nuclear force [El El El]- In the calculation, it is assumed that, at 
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a given density and temperature, matter is arranged in a periodic structure through¬ 
out a sufficiently large region of space for a unit cell to be identified, in which the 
microscopic and bulk properties of the matter are calculated. The calculation is 
performed in cubic cells with periodic boundary conditions and assuming reflection 
symmetry across the three Cartesian axes. The required reflection symmetry allows 
us to obtain solutions only in one octant of the unit cell, which reduces significantly 
the computer time. The only effect of confining ourselves to 1/8 of the cell is that 
we can only consider triaxial shapes. It is expected that the absolute minimum of 
the free energy of a cell containing A nucleons is not going to be particularly pro¬ 
nounced and there will be a host of local minima separated by relatively small energy 
differences. In order to systematically survey the “shape space” of all nuclear config¬ 
urations of interest, the quadrupole moment of the neutron density distributions has 
been parametrized, and those parameters constrained. It is expected that the proton 
distribution follows closely that of the neutrons. The minimum of the free energy in 
a cell at a given particle number density, temperature and a proton fraction is sought 
as a function of 3 free parameters, the number of particles in the cell (determining the 
cell size) and the parameters of the quadrupole moment of the neutron distribution 
/3,j. We selected four different interactions, SkM* [13], SLy4 [in], NRAPR [T7| and 
SQMC700 [iHI, based on their overall performance in modeling of a wide variety of 
nuclear matter properties PI- 

3 The RMF approximation 

We use relativistic mean field nuclear parametrizations with constant couplings and 
nonlinear terms, FSU [20] and lUFSU 125, and with density-dependent couplings, 
TW [22] . The FSU model is too soft and does not describe a 2 Mq neutron star, 
however it is expected to describe well the crust |7|. 

The equations are derived from the Lagrangian density, where we consider a sys¬ 
tem of nucleons and electrons, and is given by 



where the nucleon Lagrangian reads 


A = [iaD^ - M*] t/i. 


( 1 ) 


with 



( 2 ) 

(3) 


M* = M-r^a 
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and the electron Lagrangian is given by 


Ce = 'tpe bti - rUe] t/'e- (4) 

The meson and electromagnetic Lagrangian densities are 

^ {di^ad^a - mla^ 

= I 

r — _ 1 p 

^ p'^ 

^nl = 

where dnUJy di/COn, d^Pj^ c^j/p^ rp(p^ x Pj^) and F^i, d^Ay dyA^. 

Light clusters 

To add light clusters {d =^H, t =^H, a =^He, h =^He) to our system, we need to 
increase its degrees of freedom. The Lagrangian density becomes 

C = Ci + La Ld -\- Lg -\- Lfj + L^ + Tp + + L^i- (5) 

i=p,n,t,h 

The a particles and the deuterons are described as in [23] : 

£. = (6) 

£.i = (?) 

The effective masses of the clusters, M*,j = d,t,h,a and their binding energies 
are given in m- 

4 The Thomas-Fermi and the Coexisting Phases 
approximations 

We use the Thomas-Fermi approximation to describe the non-uniform npe matter 
inside the Wigner-Seitz unit cell, that is taken to be a sphere, a cylinder or a slab in 
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three, two, and one dimensions [2S1Z!. In this approximation, npe matter is assumed 
locally homogeneous and at each point its density is determined by the corresponding 
local Fermi momenta. In 3D we consider spherical symmetry and in 2D we assume 
axial symmetry around the z axis. The helds are assumed to vary slowly so that 
baryons can be treated as moving in locally constant fields at each point [3]. In this 
approximation, the surface effects are treated self-consistently. 

In the Coexisting Phases (CP) method, matter is organized into separated regions 
of higher and lower density, the higher being the pasta phases, and the lower a 
background nucleon gas. The interface between these regions is sharp, and it is taken 
into account by a surface term and a Coulomb one in the energy density [I2]. This 
method is much faster than the TF one, and the Gibbs equilibrium conditions are 
used to get the lowest energy state. 


5 The thermodynamical spinodal calculation 


Matter is stable with respect to fluctuations in density and composition, under condi¬ 
tions of constant volume and temperature, when the free energy density is a convex 
function of the proton and neutron densities [25] . These densities are associated with 




the chemical potentials pn = 

and Pp = 

dT 

dpp 

by 

/ 

d‘^T 

\ 

C = 

d^Pp 

dppdpn 






\ dpndpj^ 

, d'^pn 

/ 


9pp 

Ofjjp 


9pp 


\ dpn dpn / 


The eigenvalues of this matrix are given by 


A 




hTr(C) ± yD:(C)2 - 4Det(C) 


and the eigenvectors by 


5pt 


A 


_ d^Jin 

^ dp„ 


dfip 

dpn 


The thermodynamical spinodal region is then defined as the surface in {p,yp,T) 
space where the curvature of is < 0. The intersection with the Equation of State 
gives the transition density to uniform matter. 


6 Results 

We show in Table H] the nuclear matter properties for the four Skyrme models, 
NRAPR, SQMC700, SkM* and SLy4, and for the three RMF parametrizations, FSU, 
lUFSU and TW, we use in this study. 
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Model 

Po 

B/A 

K 

E 

^sym 

L 

NRAPR 

0.16 

-15.85 

226 

33 

60 

SQMC700 

0.17 

-15.49 

222 

33 

59 

SkM* 

0.16 

-15.77 

217 

30 

46 

SLy4 

0.16 

-15.97 

230 

32 

46 

ESU 

0.15 

-16.3 

230 

33 

60.5 

lUFSU 

0.155 

-16.40 

231 

31 

47 

TW 

0.15 

-16.3 

240 

33 

55 


Table 1: Nuclear matter properties at saturation density po (energy per particle B/A, 
incompressibility K, symmetry energy Egym and symmetry energy slope L) for the 
models studied. All the quantities are in MeV, except for po) given in fm“^. 


In Figure [U we show the free energy per particle as a function of the density, 
for the FSU interaction, T = 4 MeV and yp = Pp/Pn = 0.3. We show the results 
for homogeneous matter (red) and CP (blue) calculation with (dashed) and without 
(solid line) clusters. The TF results are shown with points. We can observe that F 
is lowered when pasta is present, making these states the most stable ones. The light 
clusters are only present for very small densities, and will start melting for p > 0.001 
fm“^, and their presence lowers the free energy. 

Phase transitions can be classified as first or second (or higher) order (see e.g. 
[26]). For a first-order phase transition, at least one of the first derivatives of the 
free energy with respect to one of its variables is discontinuous: S = —^\n,v,... 
P = This discontinuity produces a divergence in the higher derivatives 

like the specific heat Cy = T^\v = —or the incompressibility A'(po) = 
where Esnm is the energy per particle of symmetric nuclear matter. 
For a phase transition of second (or nth order), the first derivatives of the free energy 
are continuous; however, second (or nth order) derivatives, like the specific heat or 
the susceptibility, are discontinuous or divergent. 

In Figure [21 we plot the total pressure as a function of the density. We can observe 
that the homogeneous matter pressure curvature is partially removed with the TF and 
CP calculations, though CP gives a larger correction, which may not be so realistic. 
The intermediate transitions (between shapes) are given by a discontinuity in the total 
pressure for the TF calculation. In the 3DHFEOS approximation, a discontinuity is 
shown for the transition to homogeneous matter, and is circled in the right plot. 

In Figure [31 we plot the baryonic pressure as a function of the density. Like 
in Fig. El the homogeneous matter baryonic pressure curvature is partially removed 
within all methods, TF, CP, and 3DHFEOS, but the baryonic pressure still exhibits 
a region with a negative incompressibility. The electrons are responsible for making 
the total pressure positive. 
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Figure 1: Free energy per baryon as a function of the density for the FSU interaction, 
T = 4 MeV and Up = 0.3. Results with and without pasta, and including or not clus¬ 
ters, are shown. The effect of these aggregates are only seen for very small densities 
(inset panel). 



P p 


Figure 2: Total pressure as a function of the density for the FSU (left) and Skyrme 
(right) parametrizations, T = 4 MeV and tjp = 0.3. Results for the Thomas-Fermi 
(points). Coexisting Phases (dashed) and homogeneous matter (solid line) calcula¬ 
tions are shown in the left panel. In the right panel, the transitions to homogeneous 
matter are circled. (Right plot taken from |14j.) 
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Figure 3: Left: Baryonic (yellow) and electronic (blue) pressure components as a 
function of the density for the FSU interaction for T = 4 MeV and Up = 0.3. Results 
for the Thomas-Fermi (points), Coexisting Phases (dashed) and homogeneous matter 
(solid line) calculations are shown. Right: (Plot taken from [T3j) Baryonic pressure 
for the Skyrme interactions in the 3DHFEOS (points) and uniform matter (lines) 
calculations, for T = 2,4, 6 MeV, and i/p = 0.3. 


Let us now discuss the evolution of the geometrical pasta shapes with the density. 
In Figure m we show the neutron density distribution for the SQMC700 interaction 
and T = 2 MeV and yp = 0.3 at the onset density for each shape we found. The pasta 
configurations shown here appear for all the Skyrme models, but the threshold density 
changes somewhat. At 0.042 fm“^ a new shape, cross-rod, normally not considered in 
pasta calculations, appears. This was due to the freedom used to find the equilibrium 
configuration. 

In Figure El we show the transition densities between the pasta configurations 
for the FSU interaction (left) and the Skyrme models (right). The slabs are omitted 
(except for T = 4 MeV) for the FSU parametrization and occupy the widest density 
range in Skyrme models considered. However, the slab geometry is present in other 
RMF parametrizations, as we can see in Figure El In FSU, the difference in the free 
energy between slabs and tubes is less than 10“^ MeV. Therefore, the stable geome¬ 
tries depend on the parametrizations considered, and the properties that influence 
them should be investigated. We can also observe from the left panel that the density 
range of each shape decreases with increasing temperature. Looking at the right panel 
of Fig. El and as already seen in Fig. IH one more shape, cross-rods, is obtained with 
3DHFEOS calculations. Also, the shapes and their range of densities are almost not 
affected by the choice of the Skyrme parametrizations, and this is probably due to 
the fact that we are using a number of Skyrme models that give very close predictions 
for nuclear matter properties (as seen in Table ED . 

We have discussed so far the pasta configurations and their range of densities for a 


7 








p [fm 


0.020 

0.032(2) 

0.042(2) 

0.052(2) 

0.090(2) 0.108(2) 0.114(1) 

• 

• 

• 


• • » 

• 

• 

• 

• 

m m m 

• 

§ 

• 

■ 

■ ■ ■ 

• 

• 

• 

• 

■ ■ ■ 


Figure 4: (Figure taken from [Bj) The neutron density distribution for the SQMC700 
model, and T = 2 MeV and Hp = 0.3, is shown at the density corresponding to the 
onset of each phase, known with the uncertainty given in brackets. Blue (red) color 
indicates the bottom (top) of the density scale: 0.001 (dark blue) - 0.02475 (light 
blue) - 0.0485 (green) - 0.07225 (light orange) - 0.095 (red) fm“^. 


TF approximation 


3DHFEOS approximation 
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Figure 5: Transition densities between pasta formations. Left: FSU and T = 4 — 10 
MeV. Right (data taken from [Ej): Skyrme interactions and T = 2 MeV. 
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Figure 6: Transition densities between pasta formations for TW and T = 0 — 8 MeV 
(data taken from [23]). 

fixed proton fraction. Are there any differences if we consider /3—equilibrium matter? 
Let us observe Figure d where the transition densities are plotted as a function of 
the temperature for FSU and lUFSU and /3—equilibrium matter. We can see that 
the droplet-rod and rod-slab transition densities depend weakly on the temperature, 
and the melting T is different for all shapes. lUFSU has pasta up to T ~ 10 MeV, 
but for other parametrizations pasta melts at T ~ 4 — 6 MeV. In comparison with 
Up = 0.3 matter, /^-equilibrium matter only shows three stable configurations. FSU 
pasta melting T is lower (T ~ 6 MeV) in this case than in fixed Up matter (see Fig. E]). 

Finally, in Figs. IH1 and |9] we plot the transition densities to uniform matter. In 
Fig. El we compare the 3DHFEOS calculation with the four Skyrme interactions, with 
the CP and TF calculations for the FSU model. We can observe that the TF and 
CP calculations are almost coincident. The difference to the 3DHFEOS calculation is 
~ 0.015 fm“^ and decreases with increasing T. In the left panel of Eig. IHl we compare 
the results obtained with TE and TS for the TW model with the results obtained with 
the 3DHEEOS and TS for the NRAPR interaction. In the right panel, we demonstrate 
how to obtain the transition density to uniform matter within the thermodynamical 
spinodal calculation: the transition density is obtained from the crossing of the EoS 
at the temperature considered with the respective spinodal. We can see that the 
3DHEEOS and TE calculations give very compatible results, for the models here 
considered. When we compare with the thermodynamical spinodal calculation, we 
find that the difference of this method to the self-consistent calculations is ~ 0.015 
fm“^ for T = 0 MeV and decreases with increasing T. 
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Figure 7: Transition densities between pasta formations for FSU (left) and lUFSU 
(right, plot taken from [7j) as a function of T for /3—equilibrium matter. 


0.12 

A 

-1- 

3DHFEOS 

TF 

A 

o 

0.11 

- A 

k 

▲ 

A 

CP 

o 

0.1 


A 

A 

- 

0.09 

© 

© 

A 

A 

A 

- 

0.08 

0.07 

0.06 

SkM* 

SLy4 
■ FSU 

NRAPR 

SQMC700 


8 

1 1 

- 


4 6 8 10 

T (MeV) 

Figure 8: Transition densities to uniform matter. For more explanations, see text. 
(3DHFEOS data taken from [0].) 


10 
















CO 

'e 

dr 


0.12 

0.11 

0.1 

0.09 

0.08 

0.07 

0.06 


; 

'SDHFEO'S 

A 


TF 

♦ 

* 

A TS 

A 


♦ 


: A 

▲ 

A 

- 

; A 

♦ 

<] 

< 



A A 



A 


- 


A : 



A ; 

TW 



NRAPR 




0 2 4 6 8 10 

T (MeV) 



0 0.02 0.04 0.06 0.08 

Pn 


Figure 9: Left: Transition densities to uniform matter. Right: TS calculation. For 
more explanations, see text. (Left plot: NRAPR data taken from [H] and TW data 
taken from [24] . Right plot taken from no 


7 Conclusions 

The main conclusions of the present study are listed in the following. The effect of 
light clusters is weak and only noticeable at very low densities, but it was shown that 
taking clusters into account lowers the free energy. They will certainly affect trans¬ 
port properties. The density range of the pasta phase and the crust-core transition 
density decrease with increasing temperature, however the melting temperature of the 
different pasta phase geometries, or the inner crust, depend on the model properties. 
Our Skyrme models fall within the range of nuclear matter properties (especially the 
symmetry energy and its slope, see Table [Tj) predicted by a number of experimen¬ 
tal probes, and within this range, the pasta appearance and structure does not vary 
much. Both TF and 3DHFEOS calculations show jumps in the pressure that indicate 
a hrst order phase transition to uniform matter. Within the CP method, the inclu¬ 
sion of a non-homogeneous phase gives unrealistic results at low densities, though 
it predicts concordant transition densities to uniform matter. The stable geometries 
depend on the parametrizations used, and the properties that influence them should 
be investigated. It was also observed that /3-equihbrium and fixed proton fraction 
matter show different melting temperatures for pasta, and different number of con¬ 
figurations. All the methods considered in this study show a very good agreement 
with respect to the transition density to homogeneous matter. 
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